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Abstract 

The knee in the all-particle energy spectrum is scrutinized with a phenomenological 
model, named poly-gonato model, linking results from direct and indirect measurements. 
^ ' For this purpose, recent results from direct and indirect measurements of cosmic rays 

in the energy range from 10 GeV up to 1 EeV are examined. The energy spectra of 
, individual elements, as obtained by direct observations, are extrapolated to high energies 

using power laws and compared to all-particle spectra from air shower measurements. A 
cut-off for each element proportional to its charge Z is assumed. The model describes the 
knee in the all-particle energy spectrum as a result of subsequent cut-offs for individual 
elements, starting with the proton component at 4.5 PeV, and the second change of 
O ^ the spectral index around 0.4 EeV as due to the end of stable elements (Z = 92). The 

q . mass composition, extrapolated from direct measurements to high energies, using the 

poly-gonato model, is compatible with results from air shower experiments measuring 
the electromagnetic, muonic, and hadronic components. But it disagrees with the mass 
composition derived from X max measurements using Cerenkov and fluorescence light 
detectors. 

^ ' Key words: Cosmic rays; Energy spectrum; Knee; Mass composition; Air shower 

^ PACS: 96.40.De, 96.40.Pq, 98.70.Sa 



1 Introduction 

After almost one century of cosmic-ray physics, 
a statement made by V. Hess in his Nobel Lec- 
ture in 1936 is still true: From a consideration of 
the immense volume of newly discovered facts in 
the field of physics ... it may well appear . . . that 
the main problems were already solved .... This 
is far from being truth, as will be shown by one 
of the biggest and most important newly opened 
fields of research, . . . that of cosmic rays 0] . 
Starting with very simple instruments flown in 
manned balloon gondolas in 1912 (2| the field de- 



veloped to present-day experiments, using mod- 
ern particle physics detectors at ground level, in 
the upper atmosphere, and in space to detect 
cosmic-ray particles. Many properties of cosmic 
rays have been clarified by these measurements, 
but the origin of the relativistic particles is still 
under discussion. Since charged particles are de- 
flected by the irregular galactic magnetic fields, 
the arrival direction cannot be used to identify 
their sources, but measurements of their energy 
spectra allow to draw conclusions about their 
origin. 

The energy spectra for individual elements 



1 



are obtained directly with satellite and balloon- 
borne experiments at the top of the atmosphere 
up to several 10 13 eV and for groups of elements 
up to 10 15 eV. Due to the fast decreasing flux, 
measurements at higher energies require large 
detection areas or long exposure times, which 
presently can only be realized in ground-based 
detector systems. These experiments measure 
extensive air showers, generated by interactions 
of the high energetic cosmic rays with the nuclei 
in the atmosphere. 
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Figure 1: Cosmic-ray energy spectrum, taken from 



The primary energy spectrum extends over 
many orders of magnitude from GeV energies 
up to at least 10 20 eV. It is a steeply falling 
power law spectrum with almost no structure as 
can be inferred from Figure |l[ The only promi- 
nent features are a change in the spectral in- 
dex at about 4 PeV, generally called the knee, 
a slight steepening around 400 PeV, designated 
as 2nd knee in Figure |], and a flattening at the 
highest energies around 10 EeV, called the an- 
kle. The knee was first observed by the MSU 
group in the electromagnetic component more 
than 40 years ago Jl. Since then it has been 
confirmed by many groups in other components 
of extensive air showers, viz. the electromag- 
netic, muonic, and hadronic component, as well 



as in their Cerenkov radiation. Yet its origin is 
still under discussion and is generally believed 
to be a corner stone in understanding the origin 
of cosmic rays. 

In this article the two-knee structure is ex- 
plained as a consequence of the all-particle spec- 
trum being composed of individual spectra of 
elements with distinguished cut-offs. The first 
and prominent icnee is due to the subsequent 
cut-offs for all elements, starting with the pro- 
ton component, and the second knee marks the 
end of the stable elements (Z = 92). 

The idea of composing the all-particle spec- 
trum in the knee region out of individual elemen- 
tal spectra has been discussed in the literature, 
e.g. Webber sketches the icnee in the energy 
spectrum as result of a rigidity dependent cut- 
off for individual elements || or Stanev et al. 
calculate the all-particle spectrum as sum spec- 
trum for groups of elements from hydrogen to 
iron Q. 

In the following the phenomenological model 
mentioned is described, for which the Greek no- 
tation poly gonatof\ is chosen. In section || and 
|] an overview is given of existing direct mea- 
surements of individual nuclei up to ultra-heavy 
elements which are important to explain the sec- 
ond icnee, as well as of indirect measurements of 
the all-particle spectrum. 

Preliminary results of this work have been 
presented earlier Since then, some parame- 
ters have changed due to inclusion of recent re- 
sults and application of more detailed analysis 
procedures. 



The 
Model 



Phenomenological 



Cosmic-ray particles are most likely accelerated 
in strong shock fronts of supernova remnants 
by diffusive shock acceleration. The mechanism 
first introduced by Fermi [|| has been improved 
by many authors, see e.g. |9|, 10, 11, |l2|, In 
these models particles are deflected by moving 
magnetized plasma clouds, the particles cross 
the shock front many times, and thereby gain 
energy up to the PeV region and achieve a non- 
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thermal energy distribution. The Fermi accel- 
eration at strong shocks produces a power-law 
spectrum close to that observed. 

During the diffusive propagation of the parti- 
cles through the galaxy on a random path de- 
flected by turbulent magnetic fields, the energy 
spectrum of the particles is modified. The differ- 
ences between source spectra and observed spec- 
tra have been explained due to nuclear spalla- 
tion, leakage from the galaxy, nuclear decay, ion- 
ization losses, and, for low energies, solar mod- 
ulation, see e.g. [[□, Il4|, |i~5|. 



The poly-gonato model parametrizes the en- 
ergy spectra of individual elements, based on 
measured data, assuming power laws and taking 
into account the solar modulation at low ener- 
gies. 

At energies in the GeV region the solar mod- 
ulation of the energy spectra is described using 
the parametrization 
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adopted from J23| . N is a normalization con- 
stant, 7z the spectral index of the anticipated 
power law at high energies (see below), E = 
Eq/A the energy per nucleon, uia the mass of 
the nucleus with mass number A, and M the 
solar modulation parameter. A value M = 
750 MeV is used for the parametrization. 

Above about Z ■ 10 GeV, the modulation due 
to the magnetic field of the heliosphere is negli- 
gible and the energy spectra of cosmic-ray nu- 
clei are assumed to be described by power laws. 
The finite life time of supernova blast waves lim- 
its the maximum energy E max to be obtained 
during the acceleration process. The latter be- 
comes inefficient at an energy Ez oc Z, where 
Z is the nuclear charge of the particle, and a 
complete cut-off is expected at the maximum 
energy E max , also proportional to Z, see e.g. 
||, 0, |l|, Therefore, a cut-off or at least 
a change of the slope in the spectrum for the 
individual species at an energy Ez oc Z is ex- 
pected. 

In the leaky box model for cosmic-ray prop- 
agation particles are contained in a well de- 



fined volume of the galaxy, but have a small 
probability to escape from this volume. The 
gyromagnetic radius for a particle with charge 
Z, energy E\§ in PeV, and momentum p in a 
magnetic field B is given by R = p/eZB as 
(10 pc) E15/ZB (pG) ppfl . In a simple picture of 
the propagation low Z particles are more likely 
to escape from the galaxy as compared to parti- 
cles with high Z at the same energy due to their 
larger gyromagnetic radii. 

In diffusion models, the global toroidal mag- 
netic field in the galaxy disturbs the random 
walk of particles in the stochastic magnetic fields 
and causes a systematic drift or Hall diffusion 
plf . This drift is rapidly rising with energy. 
It dominates at the knee the slowly increasing 
usual diffusion and leads to a strong leakage 
from the galaxy. Diffusion and drift of different 
nuclei at ultrarelativistic energies depend on the 
energy per unit charge Eq/Z, and consequently 
lead to a rigidity dependent cut-off for individ- 
ual elements. 

To summarize, both scenarios — acceleration 
and propagation — entail that the cut-off en- 
ergy for each individual element depends on its 
charge Z. Inspired by these theories the fol- 
lowing ansatz is adopted to describe the energy 
dependence of the flux for particles with charge 
Z 



7c ~ Iz 
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dE 
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The absolute flux and the spectral index 7^ 
quantify the power law. The flux above the cut- 
off energy is modeled by a second and steeper 
power law. j c and e c characterize the change in 
the spectrum at the cut-off energy Ez- Both 
parameters are assumed to be identical for all 
spectra, 7 C being the hypothetical slope beyond 
the knee and e c describes the smoothness of the 
transition from the first to the second power law. 
e c = 1 corresponds to a smooth change in the 
range of about one decade of energy, larger val- 
ues describe a faster change, e.g. e c = 4 leads to 
a change within 1/5 of a decade. 

To study systematic effects, instead of a com- 
mon spectral index for all elements above the 
cut-off energy also a constant difference A7 be- 
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tween the spectral indices below and above the 
knee is tried and the spectrum for an element 
with charge Z is assumed as 
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3 Direct Measurements 

3.1 Elements with Z < 28 

Many direct measurements of the energy spectra 
of individual cosmic-ray nuclei have been per- 
formed. Figures [| to f| show compilations of 
results for protons, helium and iron nuclei. 



From an astrophysical point of view a rigid- 
ity dependent cut-off Ez oc Z is the most likely 
description, as motivated above. But, since the 
origin of the knee is still under discussion, also 
other explanations for the change of the spectral 
slope are possible. Some theories still favour a 
particle physical origin due to a change of the 
interactions in the atmosphere (see e.g. |22|j ). 
Such effects would probably depend on the en- 
ergy per nucleon Eq/A and lead to a nuclear 
mass dependent cut-off Ez oc A. The last possi- 
bility considered in this article is, that the bend 
in the individual spectra appears at a constant 
energy. In the poly-gonato model all three as- 
sumptions to parametrize the cut-off energy are 
scrutinized. 



E p ■ Z ; rigidity dependent 
Ez = { E p ■ A ; mass dependent 
E p ; constant 



(4) 



It will be investigated with measured data, 
which ansatz describes the data best. In the fol- 
lowing, only "poly-gonato model" refers to the 
rigidity dependent hypothesis. 

Knowing the flux d&z /dEo(Eo) for the 
species with charge Z, the flux of the all-particle 
spectrum is obtained by summation over all 
cosmic-ray elements 



dK 
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This equation contains the parameters and 
72 for each element and the common parame- 
ters, E p , e c , and j c or A7. and jz are ob- 
tained from direct measurements of individual 
nuclei. The remaining three parameters are de- 
rived from a fit to the all-particle spectrum as 
obtained by indirect measurements. 
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Figure 2: Differential energy spectrum for protons. 
The best fit to the spectrum according to a power 
law is represented by the solid line, the bend (dot- 
ted line) is obtained from a fit to the all-particle 
spectrum, see section |^. The all-particle spectrum 
is shown as dashed line for reference. Results are 
shown from AMS ||, BESS ||, CAPRICE |§, 
HEAT [§7|, Ichimura et al. |§, IMAX §§, JACEE 
[30j, Kawamura et al. 0, MASS pi, Ormes et al. 
[33 1, Papini et al. Jl, RUNJO B J3j| , Ryan et al. 
[36|, Smith et al. §7], SOKOL |3jf, Webber et al. 
[39 1, and Zatsepin et al. [H. 



The differential flux for the most abundant 
cosmic-ray element — hydrogen — multiplied 
by Eq- 5 is plotted in Figure |2| versus the par- 
ticle energy Eq. The data show a decreasing 
flux as function of energy. The straight line rep- 
resents the fit according to a power law. The 
errors as specified by the different experiments 
have been taken into account. Best fit values are 
$° = (8.73 ± 0.07) • I0~ 2 (m^srTeV)" 1 and 
7 P = -2.71 ±0.02 (solid curve) with x 2 /d-o.f. = 
3.4. The errors specify the statistical uncertain- 
ties. The cut-off shown in the figure as dot- 
ted line will be discussed in section |5| when fit- 
ting the all-particle spectrum within the poly- 
gonato model. The all-particle spectrum ob- 
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tained is shown for reference as well (dashed 
curve) . 

Recent measurements in the 100 GeV region, 
i.e. AMS, BESS, CAPRICE, and IMAX, obtain 
lower flux values as have been reported by earlier 
experiments. These values had not been avail- 
able for an earlier compilation by Wiebel-Sooth 
et al. |l(|. Therefore, the flux in this reference 
is about 30% larger than the present result and 
the spectrum is slightly steeper. 

When taking the experimental data in Fig- 
ure |2] at face value the reader might imagine a 
change of slope around 10 4 GeV. However, the 
new measurements do not confirm such former 
speculations. Also indirect measurements give 
no hint for deviations from a power law. De- 
riving the primary proton spectrum via the in- 
vestigation of unaccompanied hadrons at moun- 
tain level, the EAS-Top experiment finds no ev- 
idence for a structure in the spectrum up to 
10 5 GeV El. 
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Figure 3: Differential energy spectrum for helium 
nuclei. The best fit to the spectrum according to a 
power law is represented by the solid line, the bend 
(dotted line) is obtained from a fit to the all-particle 
spectrum, see section g. The all-particle spectrum is 
shown as dashed line for reference. Used are values 
from Anand et al. @, BESS @, CAPRICE @, 
HEAT [§7), Ichimura et al. @, IMAX |§, JACEE 
[30], Kawamura et aL§|], MASS_jf|, Ormes et al. 
Ill, Papini et al. fM||7RIC H [o TRUNJOB ' 
Ryan et al. pi Smith et al. M, SOKOL & and 
Webber et al 

Figure [3| shows the differential energy spec- 
trum for helium nuclei as a function of the to- 



tal particle energy. The parameters for the 
spectrum, obtained by a least-square fit, are 
<$> He = (5.71 ± 0.09) • 10" 2 (m 2 ssrTeV) _1 and 
-/ He = -2.64 ± 0.02 with x 2 /d-o.f. = 3.5. The 
spectral index agrees with the value obtained by 
Wiebel-Sooth et al., but the absolute flux is now 
about 25% lower due to contributions of recent 
measurements, both at the low and high energy 
end of the spectrum. 
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Figure 4: Differential energy spectrum for iron nu- 
clei. The best fit to the spectrum is represented 
by a solid line, the all-particle spectrum is shown 
as dashed line for reference. Plotted are measure- 
ments according to CRN H], Engelmann et al. p5|, 
Hareyama et al. [Q, Ichimura et al. p8| , JACEE 
p7[ , Juliusson et al. Q, M inagawa et al. ||49|, Orth 
et al. IF" 
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As an example for heavy elements the en- 
ergy spectrum for the most prominent repre- 
sentative — iron nuclei — is presented in Fig- 
ure |j. Results from detectors which are able to 
resolve individual elements are given together 
with experiments, where the charge resolution 
allows only to specify the flux for the iron 
group, e.g. JACEE, RUNJOB, Simon et ah, 
and SOKOL. The best fit to all measurements 
shown above 1 TeV yields for the power law 
$° e = (2.04 ±0.26) • 10~ 2 (rn^ssrTeV)- 1 and 
7 Fe = -2.59 ± 0.06 with x 2 /d.o.f. = 0.9. If the 
result is compared to the earlier compilation, 
the present absolute flux is about 15% larger, 
while the spectral index agrees within the spec- 
ified errors with the old value. The suppression 
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Figure 5: Spectral index 7^ versus nuclear charge 
see Table ^. The solid line represents a three 
parameter fit according to Equation (^), the dashed 
graph a linear fit. 



of the flux at energies below 1 TeV is due to 
the modulation in the heliosphere. It can be in- 
ferred from the figure, that the measurements 
are parametrized well using Equation (|l]). 

The fit parameters for protons, helium, and 
iron nuclei are listed together with results from 
Wiebel-Sooth et al. for the elements hydrogen 
to nickel in Table [7] in the appendix. The spec- 
tral indices for these elements are plotted in Fig- 
ure [5] versus the nuclear charge. The data show a 
trend towards smaller values with increasing Z. 
Such a Z-dependence of 7^ could be explained 
by charge dependent effects in the acceleration 
or propagation process. 

Theories using a nonlinear model of Fermi accel- 
eration in supernovae remnants predict a more 
efficient acceleration for elements with a large 
mass to charge ratio as compared to elements 
with a smaller A/Z ratio, see e.g. [53]. Conse- 



quently, it is expected that elements with higher 
A/Z have a harder spectrum at the source. 
The energy spectra observed at earth are modi- 
fied during propagation of the particles through 
the galaxy. Some authors include reacceleration 
by weak interstellar shocks in the standard leaky 
box model, e.g. [05, 2C, pM. Like the primary ac- 



celeration also the reacceleration could be more 
efficient for high-Z-nuclei. 

To estimate the fluxes of these ultra-heavy 



elements at high energies the parametrization 

- 7Z = A + B ■ Z c (6) 

is used to describe the Z dependence of the spec- 
tral indices and to extrapolate them to higher 
values. To study systematic effects of the ex- 
trapolation two approaches are used, a linear 
function (C = 1) and a non-linear extrapola- 
tion, using all three parameters. 

The dashed line in Figure || represents the 
best fit of a linear parametrization, exhibiting 
a decreasing spectral index as function of the 
nuclear charge. The data shown in the figure 
exhibit some curvature which suggests to in- 
troduce the additional degree of freedom. If 
the parameter C in Equation (0) is used as 
free parameter, the solid line in Figure [| is ob- 
tained. The parameters for both trials are listed 
in Table |], both fits result in about the same 
X 2 /d.o.f. ~ 2.1. The values for the non-linear 
approach will be corroborated below by an in- 
dependent fit to the all-particle spectrum. 

3.2 Ultra-heavy elements 
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Figure 6: Relative abundance of cosmic-ray ele- 
ments (Z > 28) normalized to Fe= 1 from ARIEL 6 
||, Fowler et al. [||, HEAO 3 |7), SKYLAB §§, 
TREK/MIR fH, Tueller et al. and Israel et al. 
[ jSP] , as well as UHCRE JsT|] . The energy is around 
1 GeV/nucleon. 

For ultra-heavy elements {Z > 28) data ex- 
ist only at relative low energies around a few 
GeV/nucleon as already mentioned. Figure |6] 
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Table 1: Parameters of Equation (g) for the linear and non-linear extrapolation of jz. 



ABC 
linear 2.69 ±0.12 (-2.07 ± 1.05) • 1(T 3 =1 

non-linear 2.70 ± 0.19 (-8.34 ± 4.67) • 10~ 4 1.51 ± 0.13 



shows a compilation of the relative abundance 
from copper [Z = 29) up to uranium (Z = 92), 
as measured by several experiments on space 
crafts and balloons. The data are normal- 
ized to Fe= 1, the threshold is about 0.5 to 
1 GeV/nucleon. Some authors give only results 
for groups of elements, this is indicated by hor- 
izontal error bars. 

The experiments ARIEL 6 §5|], HEAO 3 [El, 
as well as Tueller et al. and Israel et al. [ pOfl 
quote abundances relative to iron. Only relative 
abundances for elements Z > 70 are reported by 
Fowler et al. §§, SKYLAB [EJ, TREK/MIR 



H, and UHCRE Q. The results of the latter 
group have been normalized to ARIEL 6. This 
detector could resolve individual elements up to 
Z = 48, and even charged nuclei above. The 
range 70 < Z < 80 has been used to match 
the abundances for Fowler et al., SKYLAB, and 
UHCRE. For TREK the interval 72 < Z < 80 
has been utilized. 

The results of all experiments show about the 
same structure for the relative abundances. For 
elements with Z > 80 deviations are visible. 
Due to the very low flux only a few (< 10) nuclei 
have been detected during a typical mission and 
the experiments reach their limit for statistically 
reliable results. 

The relative abundances shown in Figure ^ 
are the only data available for ultra-heavy ele- 
ments. No data about individual spectra are at 
disposal. Hence, the spectral indices for these 
elements have to be estimated and are obtained 
by extrapolating parametrization (^) to large Z. 
The mean relative abundances for the elements 
are calculated from the data given in Figure ^, 
converted to absolute flux values and extrapo- 
lated to higher energies in the following manner. 
The solar modulation is taken into account using 
Equation (||) and the fluxes are extrapolated to 
1 TeV/nucleus assuming at high energies power 
laws with the extrapolated spectral indices. 
For Z > 48 the detectors could only resolve even 



charged elements, i.e. the values shown are the 
sum of two elements. For the extrapolation they 
are split between odd and even atomic numbers 
using the ratio 1:5. This ratio is a typical value 
for the solar system composition, see Figure |7]. 
The values thus obtained for <& z and jz, using 
the non-linear extrapolation to obtain 7^, are 
listed in Table 0. 

3.3 Comparison with solar system 
composition 
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Figure 7: Absolute flux &° z 01 cosmic-ray elements 
at -Eq = 1 TeV/nucleus versus nuclear charge ac- 
cording to Table 0. Two versions for the extrapola- 
tion of "fz for ultra-heavy nuclei are indicated, see 
text. For comparison, solar system abundances are 
shown ||6^| . The relative values for the solar system 
have been normalized to the absolute flux of silicon 
[Z = 14) in cosmic rays. The ordinate accounts for 
cosmic rays only. 

The absolute flux of all cosmic-ray elements 
at 1 TeV/nucleus, as listed in Table ^, is plot- 
ted in Figure |?] versus the nuclear charge. For 
the ultra-heavy elements the linear extrapola- 
tion (C = 1) yields slightly smaller flux val- 
ues as compared to the non-linear description 
(C > 1). For comparison, solar system abun- 
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dances are shown as well. The relative abun- 
dances for the solar system are normalized to 
the absolute flux of silicon in cosmic rays. The 
general structure of both distributions is about 
the same and the differences between even and 
odd nuclear charge numbers are mostly visible 
in both compositions. However, there are differ- 
ences which will be discussed next. 

For cosmic rays the low abundance "valleys" 
in the solar system composition around Z=4, 
21, 46, and 70 are not present. This is usually 
believed to be the result of spallation of heav- 
ier nuclei during their propagation through the 
galaxy. Hydrogen, helium, and the CNO-group 
are suppressed in cosmic rays. This has been 
explained by the high first ionization potential 
of these atoms |3jl] or by the high volatility of 
these elements which do not condense on inter- 



stellar grains [S4|. Which property is the right 
descriptor of cosmic-ray abundances has proved 
elusive, however, the volatility seems to become 
the more accepted solution J65[ . 

The iron group and the ultra-heavy elements 
are more pronounced in cosmic rays as compared 
to the solar system. Especially the r-process 
elements beyond xenon (.2=54) are enhanced, 
partly due to spallation products of the plat- 
inum and lead nuclei (Z=78, 82). For the lat- 
ter direct measurements at low energies around 
1 GeV/n yield about a factor two more abun- 
dance as compared to the solar system and a 
factor of four for the actinides thorium and ura- 
nium (Z=90, 92) j66|. This has been attributed 
to the hypothesis that cosmic rays are accel- 
erated out of supernova ejecta-enriched matter 

& 

The differences for ultra-heavy elements be- 
tween cosmic rays and solar system abun- 
dances seen in Figure ^ are much bigger as 
compared to the measurements at energies of 
1 GeV/nucleon. At low energies the flux is 
strongly suppressed due to the solar modula- 
tion, while at 1 TeV/particle the effects of the 
heliospheric magnetic fields are almost negligi- 
ble. Together with the enhancements discussed 
before, this accounts for the differences seen in 
Figure 0. The overall distribution from hydro- 
gen to uranium is much natter for cosmic rays. 
While the solar system abundances cover a range 



of 11 decades, the cosmic-ray distribution ex- 
tends over 6 decades in flux only. 

4 Indirect Measurements 

Many groups published results on the all- 
particle energy spectrum from indirect mea- 
surements. Several experiments detect the 
main components of extensive air showers 
most commonly the electromagnetic component, 
but also muons and hadrons are investigated. 
The Cerenkov photons produced by relativis- 
tic shower particles and the fluorescence light of 
nitrogen molecules in the atmosphere induced 
by air showers are also utilized. The detectors 
are located at various atmospheric depths cor- 
responding to altitudes from 50 m up to 4370 m 
a.s.L. The experiments used in this compila- 
tion and their respective measuring techniques 
as well as the atmospheric overburden are listed 
in Table g. 

Some of the experiments present different re- 
sults, based on different hadronic interaction 
models used to interpret the data. Detailed 
studies of air shower properties, especially con- 
cerning the hadronic component, have been per- 



formed by the KASCADE group §g Com- 
parison to simulations, using the program COR- 



SIKA [71 1 with several high-energy hadronic 
interaction models implemented revealed, that 
the QGSJET model presently is the best model 
to describe air shower data. Other groups ob- 
tained similar results, see e.g. |7|,0]. Hence, if 
more than one interpretation of measurements 
is given, the energy spectrum derived with the 
QGSJET model is used. 

The results of the groups are presented in Fig- 
ure ||. The overall agreement between the exper- 
iments is quite good, the differential fluxes mul- 
tiplied by -E 2,5 agree within a factor of two. All 
experiments exhibit a similar shape of the spec- 
trum despite their different absolute normaliza- 
tion. This is remarkable since different compo- 
nents of air showers are investigated by various 
groups at different atmospheric depths, i.e. dif- 
ferent types of particle interactions in the atmo- 
sphere are probed and the showers are sampled 
at different stages of their longitudinal develop- 
ment. 
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Figure 8: All-particle energy spectra obtained from direct and indirect measurements, for references see Table | 
and text. The sum spectra for individual elements according to the poly-gonato model are represented by the 
dotted line for 1 < Z < 28 and by the solid line for 1 < Z < 92. Above 10 8 GeV the dashed line gives the 
normalized average spectrum. 



The all-particle spectrum obtained by direct 



measurements from JACEE g7|], RUNJOB |pj}|, 
SOKOL [38], and Grigorov et al. p§ agrees 



with the indirect measurements in the region of 
overlap, i.e. from below 100 TeV to 1 PeV par- 
ticle energy. 
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Figure 9: All-particle energy spectrum at very high 
energies according to the Fly's Eye experiment. The 
dashed line represents a fit by Bird et al. jp8|. 



The knee at about 4 PeV is clearly recogniz- 
able in the spectrum. Assuming this bend be- 
ing caused by the cut-off of the proton com- 
ponent, the galactic component extends in the 
poly-gonato model up to Zjj ■ 4 PeV « 0.4 EeV. 
The energy spectrum multiplied by Eq is shown 
in Figure ^ for high energies as reported by 
the Fly's Eye group. Indeed, a change in the 
spectral slope around 4 • 10 8 GeV is visible, 
the dashed line represents a fit taken from the 
Fly's Eye publication [88]. A similar structure 
has been observed by the experiments AKENO 
and Haverah Park, the spectral indices obtained 
around 400 PeV by the three experiments are 
summarized in Table ||. The change coincides 
well with the anticipated cut-off of the heaviest 
nuclei of the galactic component. 



5 Cosmic— ray energy spectrum 

The different absolute normalizations between 
the experiments, apparent in Figure ||, are prob- 
ably caused by uncertainties in the energy cal- 
ibration. The latter are quoted by the experi- 
ments to be in the order of 10% to 20%. That 
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Table 2: Air shower experiments and components measured to derive the primary energy spectrum: e: elec- 
tromagnetic, /i: muonic, h: hadronic, C: Cerenkov, and F: fluorescence light. The particle thresholds are given 
for the muonic and hadronic components. In addition, the atmospheric overburden [g/cm 2 ] and the shift of the 
energy scale are listed. 



Experiment 



C F g/cm Energy shift 



AKENO (low energy) 
BLANCA @ 
CASA-MIA WEH 



DICE |7§ 
EAS-Top 
HEGRA f7| 

K AS CADE (electrons/muons 
KASCADE (hadrons/muons) |8 
KASCADE (neural network) [fTl 



MSU pi 



Mt. Norikura §|] 



Tibet §4] 
Tunka-13 @ 
Yakutsk (low energy) 



x 
x 
x 
x 
x 
x 
x 

x 
x 
x 
x 



1 GeV 

800 MeV 
800 MeV 
1 GeV 

230 MeV 
230 MeV 
230 MeV 



50 GeV 



x 
x 



930 
870 
870 
860 
820 
790 
1022 
1022 
1022 
1020 
735 
606 
680 
1020 



-47o 
4% 
4% 
1% 
-11% 
-10% 
-7% 
-1% 
-8% 
-5% 
9% 
-10% 
0% 



AKENO (high energy) |8 
Fly's Eye @ 
Gauhati [ |SS| ] 
Haverah Park @ 
HiRes-MIA ^1 



Yakutsk (high energy) [92] 



800 MeV 
1 GeV 



930 
860 
1025 
1018 
860 
1020 



-16% 
-3% 
-5% 
-10% 
-5% 
-20% 



Table 3: Spectral index around 400 PeV measured by three experiments, from p8| . 



Experiment 


energ 


y range 


spectral index 


enerj 


IY range 


spectral index 


AKENO 


10 15.7 


- 10 17 - 8 : 


3.02 ± 0.03 


10 17.8 


- 10 18 - 8 : 


3.16 ±0.08 


Haverah Park 


10 17.48 


- 10 17 - 6 : 


3.01 ±0.02 


10 17.6 


- 10 18 - 6 : 


3.24 ±0.07 


Fly's Eye 


10 17.3 


- 10 17 - 6 : 


3.01 ±0.06 


10 17.6 


- 10 18 - 5 : 


3.27 ±0.02 



means the energy scales may be shifted by this 
amount with respect to each other. This has 
been done in order to elaborate a more consis- 
tent all-experiment spectrum. 

5.1 Renormalization 

Data from experiments which start to measure 
below 1 PeV have been renormalized to fit the 



all-particle flux at 1 PeV extrapolated from the 
direct measurements. The solid line in Figure |8] 
below 1 PeV represents the average flux of these 
measurements. The renormalization factors ob- 
tained are listed in Table || for the detector sta- 
tions AKENO (low energy) to Yakutsk (low en- 
ergy). Only small shifts (< 11%) are necessary, 
all within the mentioned energy uncertainties. 
The mean value amounts to(— 2.9±1.7)%. 
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Figure 10: Normalized all-particle energy spectra for individual experiments. The renormalization values for 
the energy scale and references are given in Table || The sum spectra for individual elements according to the 
poly-gonato model are represented by the dotted line for 1 < Z < 28 and by the solid line for 1 < Z < 92. 
Above 10 8 GeV the dashed line reflects the average spectrum. 



For the experiments measuring at very high 
energies (AKENO (high energy) to Yakutsk 
(high energy) in Table ^) the energy scale has 
been checked subsequently in the following way. 
The spectra have been arranged according to 
the lower limit of the given energy range, i.e. 
HiRes, AKENO, Fly's Eye, Haverah Park, and 
Yakutsk, and then normalized to the average 
flux of all previous experiments at their lowest 
energies. The shape of the Gauhati energy spec- 
trum is systematically different from all others. 
Therefore, it has been normalized to fit the aver- 
age flux obtained by all other experiments above 
10 8 GeV. The energy shifts necessary are listed 
m Table |, they are all negative with an average 
value (-9.8 ± 2. 



Due to differences in the shape of the spectra 
of individual experiments with respect to the av- 
erage spectrum it is not always obvious how to 
normalize. This yields errors of about ±2% for 
the renormalization of the energy scales given. 
The average energy corrections in both groups 
are negative and average to (— 5.0 ± 1.6)%. This 
indicates that the experiments tend to overesti- 
mate the energy, which could be explained by 



a systematic effect in the hadronic interaction 
models applied. 

The normalized spectra are shown in Fig- 
ure 10. The spread between the different results 



has shrunk, all fluxes are contained in a narrow 
band. The similar shape becomes obvious in this 
representation, an interesting observation, since 
the energy scaling does not change the shape of 
the spectra. 

The weighted average of the normalized all- 
particle spectra has been calculated, taking the 
errors of the individual measurements into ac- 



count, the result is presented in Figure 11 and 
Table §. The errors reflect the r.m.s. values of 
the data. The shape of this multi-experiment 
all-particle energy spectrum, already seen in 
Figure [To|, becomes clearly apparent in this plot. 



5.2 Energy dependence for the cut- 
off 

The all-particle spectrum is used to determine 
the three parameters of the poly-gonato model. 
A least-square fit uses the values <&° z and 7^, 
as given in Table [?], and assumes a rigidity de- 
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Figure 11: Average all-particle energy spectrum. The line through the data represents a fit of the sum 
spectrum for elements with 1 < Z < 92 according to the poly-gonato model with rigidity dependent cut-off for 
a) common j c and b) common A7. The dotted line shows the spectrum for 1 < Z < 28. In addition, energy 
spectra for groups of elements are shown. Above 10 8 GeV the dashed line reflects the average spectrum. 

Table 4: Parameters to fit the all-particle energy spectrum for the different approaches to describe the cut-off 
behaviour, see Equations (0) to (0), with the poly-gonato model. 



cut-off: 


rigidity 


mass 


constant 






dependent 


dependent 






E z = 


E p ■ Z 


Ep-A 


E p 




E p [PeV] = 


4.51 ±0.52 


3.66 ± 0.41 


3.50 ± 0.38 


common "y c 


7c = 


-4.68 ±0.23 


-7.82 ± 1.09 


-3.06 ±0.02 




e c = 


1.87 ±0.18 


2.30 ±0.23 


1.94 ±0.51 




X 2 /d.o.f = 


0.116 


0.290 


0.086 




E p [PeV] = 


4.49 ± 0.51 


3.81 ± 0.43 


3.68 ± 0.39 


common A7 


A7 = 


2.10 ±0.24 


5.70 ± 1.23 


0.44 ± 0.02 




e c = 


1.90 ±0.19 


2.32 ±0.22 


1.84 ±0.45 




rVd.o.f. = 


0.113 


0.292 


0.088 





pendent cut-off. The two versions to model the 
cut-off behaviour are investigated, viz. a com- 
mon spectral index 7 C above the knee for all el- 
ements — Equation (Q) — or a common differ- 
ence in spectral slope A7 — Equation (||). The 
best fits obtained are shown in Figure O as solid 



lines for a common -j c (Figure 11a) and for a 
common A7 (Figure |lT|b). The parameters ob- 
tained are listed in the first column of Table ||. 

As seen in the figure, both approaches de- 
scribe the data quite well and yield low x 2 /d.o.f. 
values which are given in the table. The r.m.s. 
values of the data have been used as errors for 



the fit procedure, thus the x 2 values obtained 
are < 1. But their relative numbers are still a 
good quantity to characterize the different fit re- 
sults. The description of the cut-off behaviour 
according to Equations @ or (|3|) is not crucial 
for the results obtained. Both possibilities yield 
almost the same shape of the energy spectra as 
can be seen in Figure [ll|. To check systematic ef- 
fects, different bin sizes for the all-particle spec- 
trum have been tried and the upper limit of the 
fit range has been varied from 20 to 150 PeV. 
All fit results agree within the quoted errors. 

To derive the energy spectra for ultra-heavy 
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elements, the non-linear extrapolation for 7^ 
with C > 1 has been used. Applying the linear 
extrapolation results in a sum spectrum which 
nearly coincides with the dotted line in Figure [H] 
representing the sum spectra for 1 < Z < 28. 
The all-particle spectrum below 10 8 GeV is 
not effected by the choice for the extrapola- 
tion. Above this energy the contribution of the 
sum spectrum to the all-particle spectrum for 
1 < Z < 92 is larger for the non-linear extrap- 
olation. For example, at 5 • 10 8 GeV the sum 
spectrum contributes with about 50% for the 
non-linear and with 12% for the linear extrapo- 
lation. 

A consistency check concerning the "fz extrap- 
olation has been performed in the following way: 
The sum spectrum for 1 < Z < 92 has been 
fitted to the all-particle spectrum and the spec- 
tral indices for Z > 28 have been obtained from 
this fit, i.e. the parameters E p , e c , j c or A7 in 
Equations @ or (|3|) and the parameters A,B, 
and C in Equation (||) have been determined 
simultaneously in a six-parameter fit. The val- 
ues for both methods (7 C and A7) are consistent 
with the results listed in Table ||. The parame- 
ters for 7^ are almost identical to the values for 
the non-linear extrapolation given in Table |]. 
For this six-parameter fit the upper fit bound- 
ary has been varied from 30 PeV to 500 PeV. It 
is worth mentioning that the fit can not lower 
the exponents of the ultra-heavy elements even 
more in order to fill up the missing flux to the 
all-particle spectrum in the energy region above 
600 PeV. The boundary of the all-particle spec- 
trum at lower energies is the limiting factor. 
That means, the values for 7^ have been derived 
in two independent ways, a fit of 7^ versus Z in 
Figure ||| and in a fit of the all-particle spectrum 
with the poly-gonato model. Both fits yield the 
same results, an intriguing observation, favour- 
ing the non-linear extrapolation of -fz- 

All parameters of Equation (||) and (J3j) being 
determined, the spectra of individual elements 
can be calculated. In order to avoid confusion 
not all individual spectra, but only graphs for 
six groups of elements are plotted in Figure [H] 
for demonstration. One observes that around 
50 PeV about 50% of the particles should be- 
long to the iron group and hydrogen contributes 



only with a few % to the all-particle flux. In- 
deed, the AKENO experiment |m| found when 
determining the mass composition from electron 
and muon sizes above 20 PeV that the compo- 
sition becomes heavier with increasing energy. 
In an analysis with three mass groups (proton- 
helium, CNO, and iron) the amount of the light 
group turned out to be small, viz. in the 10% 
region or below around 50 PeV. 

The average spectrum shows some slight 
structure above 4.5 PeV, i.e. it is not a pure 
power law. In the phenomenological approach 
such a shape is attributed to the sum of 92 spec- 
tra with individual cut-off energies. 

Some authors, see e.g. [ p5]j , propose an addi- 
tional component being necessary in addition to 
the known fluxes from direct observations in or- 
der to describe the all-particle energy spectrum 
in the knee region. It may be pointed out, that 
the sum spectrum of all individual elements is 
sufficient to compose the all-particle spectrum 
up to 100 PeV, as can be inferred from Fig- 
No additional cosmic-ray com- 



10 and 11 



ures 

ponent is required in the knee region to explain 
the observed spectrum. 

In order to check the two other assumptions 
about the cut-off behaviour of the individual el- 
ement spectra — see Equation (||) — the average 
energy spectrum has been fitted using these two 
hypotheses as well. 

The results for the mass dependent cut-off us- 
ing Equations (§) and (||) are shown in Figure |l^ 
and the fit values obtained are listed in the sec- 
ond column of Table ||. Again, the resulting 
spectra are essentially independent of the ansatz 
used to describe the spectral indices above the 
cut-off. The upper end of the fit range has been 
varied from 50 to 200 PeV, yielding consistent 
results. 

The sum spectra obtained describe the over- 
all shape of the spectrum less well as compared 
to the rigidity dependent assumption. They are 
below the experimental data around 3 PeV and 
above them around 50 PeV, this is reflected in 
a worse x 2 /d.o.f. value of about a factor of 2.5 
when compared to the rigidity dependent choice. 
Also, the steepening around the individual knees 
turns out to be quite different when compar- 
ing the spectra of the elemental groups in Fig- 



13 




Figure 12: Average all-particle energy spectrum. The line through the data represents a fit of the sum spectrum 
for elements with 1 < Z < 92 according to the poly-gonato model with mass dependent cut-off for a) common 
7 C and b) common A7. The dotted line shows the spectrum for 1 < Z < 28. In addition, energy spectra for 
groups of elements are shown. 



ures |Tl] and 12. For the rigidity dependent cut- 
off the fit yields a moderate steepening, viz. a 
change of the slope of A7 = 2.10 or an index of 
7 C = —4.7 and a smooth transition between the 
power laws, ranging across slightly less than one 
decade, with e = 1.9. 

The mass dependent ansatz results in a strong 
steepening of the spectra to 7 C = —7.8 above 
the knee or a change of the spectral index of 
A7 = 5.7 within a region of about half a decade 
(e = 2.3). Taking the propagation of cosmic 
rays through the galaxy into consideration, such 
a sharp cut-off would be hard to explain on as- 
trophysical grounds. Maybe a nearby cosmic- 
ray source, as discussed in |K| , or a new type of 
interaction in the atmosphere could lead to such 
a sharp cut-off. 

Using the linear extrapolation for jz does 
not improve the x 2_vaiUe ) the spectrum around 
3 PeV can not be described satisfactorily by 
mass-dependent cut-off energies. Neverthe- 
less, a recent analysis, combining electromag- 
netic shower size spectra of several experiments 
favours a knee for individual elements propor- 
tional to their mass | |97| . 

When probing the third hypothesis for the 
cut-off energy, the constant cut-off, the up- 
per fit boundary has been changed from 20 to 
560 PeV resulting in consistent parameters, too. 
The obtained spectra are plotted in Figure ^ 



for the two assumptions on the spectral indices 
above the knee and the corresponding parame- 
ters are listed in the last column of Table For 
the common 7 C the differences between the lin- 
ear and the non-linear extrapolation for 72 are 
negligible. For a common A7 deviations above 
200 PeV are significant, as can be inferred from 
Figure p^b. The non-linear extrapolation (solid 
line) results in larger flux values as compared to 
the linear extrapolation which coincides nearly 
with the sum spectrum for 1 < Z < 28 (dotted 
line). 

The shape of the average energy spectrum is well 
described up to 200 PeV, as indicated by the low 
X 2 /d.o.f. values. But the change in the spec- 
trum at 400 PeV, seen in Figures |9| and |l4|, can 
not be described. The all-particle spectrum lies 
significantly above the measured values. A very 
smooth transition is obtained, the spectral in- 
dices change only little by A7 = 0.4 to a value 
of 7 C = —3.1 above the .knee. The width of the 
transition region, described by e = 1.9 is compa- 
rable to the one for the rigidity dependent cut- 
off. However, the ansatz of a constant cut-off 
will be disfavoured below, when discussing the 
resulting mean logarithmic mass. 

Combining all arguments discussed, the rigid- 
ity dependent cut-off is the best choice to de- 
scribe the all-particle spectrum. In the follow- 
ing the non-linear extrapolation for 72 will be 
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Figure 13: Average all-particle energy spectrum. The line through the data represents a fit of the sum spectrum 
for elements with 1 < Z < 92 according to the poly-gonato model with constant cut-off energy for a) common 
7 C and b) common A7. In addition, energy spectra for groups of elements are shown. 



used together with the hypothesis of a common 



A7 as shown in Figure |ll|b. 

The resulting spectra with a rigidity depen- 
dent cut-off for protons, He, and Fe are rep- 
resented in Figures |2] to |I| as solid lines with 
a dotted continuation in the knee region. The 
all-particle flux summed up from all individual 
elements is shown as dashed line for reference. 
In Figures |8| and |l0| the sum spectrum of indi- 
vidual elements is given for 1 < Z < 92 as solid 
line and as dotted line for 1 < Z < 28. 

The normalized spectra for experiments in- 
vestigating very high energies are shown in 



Figure 14 with the flux being multiplied by 
.Eg- The slight steepening in the spectral slope 
around 0.4 EeV can be recognized in the re- 
sults of AKENO, Fly's Eye, Haverah Park, and 
HiRes/MIA. A similar steepening occurs in the 
Yakutsk data, but at about 4 times higher ener- 
gies. Due to this unexplained systematic effect, 
the Yakutsk data are not taken into account 
when calculating the average flux. The latter 
is represented by a dashed line up to 3 EeV. 
Sum spectra for two groups of elements using a 
rigidity dependent cut-off with a common A7, 



as presented in Figure lib, are shown and it 



can be inferred, that the steepening of the spec- 
trum coincides with the cut-off of the heaviest 
nuclei according to the poly-gonato model. The 




energy E v 
reference. 



92 • E„ = 414 PeV is marked for 



10 

Energy E [ GeV ] 



Figure 14: Normalized all-particle energy spectra 
multiplied by Eq for several experiments. For nor- 
malization factors and references see Table pj. The 
dashed line represents the average flux. The sum 
spectrum according to the poly-gonato model is rep- 
resented by the solid line for all elements (.2=1-92), 
as well as for two elemental groups. The cut-off en- 
ergy Eif — 414 PeV is marked by the pointed vertical 
line. 



This may lead to the conjecture, that the sec- 
ond knee in the all-particle spectrum is caused 
by the end of the galactic component of cosmic- 
rays, i.e. by the cut-off of the heaviest elements. 
At higher energies an additional component is 
required to account for the observed flux. For 
example, at 400 PeV the sum for 1 < Z < 92 



15 



^ 10" R 



> 

o 



v 1 °' 

i 

ill 



LU 



LU 

T3 

o 

"D 



pio — 



10 — 



1 — 



10 — 



10 — 



10 — 




+ Akeno 

O BLANCA 

n KASCADE (e/m) 

$ Fly's Eye 

❖ CASA-MIA 

* KASCADE (nn) 
©HiRes/MIA 

DICE 
OMt. Norikura 



10 10 10 

Energy E [GeV] 

•w-Haverah Park 

® Tibet ASy 

A KASCADE {him) 

-J-Gauhati 

X EAS-Top 

0MSU 

>K Yakutsk 

VTunka-13 

□ HEGRA 



Figure 15: Normalized all-particle energy spec- 
tra for individual experiments compared to the all- 
particle spectrum of the poly-gonato model. The in- 
dividual results are shifted in steps of half a decade 
in flux in order to reduce overlap. 



accounts for about 55% of all particles. The lin- 
ear extrapolation for 'jz yields a spectrum which 
nearly coincides with the spectrum shown for 
1 < Z < 28. Figure 14 shows, that ultra-heavy 



all-particle spectrum in the region from 0.1 to 
1 EeV and to explain the second knee. 

In the figures presented so far it was hard 
to recognize the shape of the individual spectra 
with respect to each other. To elucidate the sit- 
uation, the normalized results of the individual 



experiments are compared in Figure 15 to the 
all-particle spectrum. The solid lines represent 
the sum spectrum up to 10 8 GeV, the dashed 
lines above represent the average spectrum. In 
order to reduce overlap between the graphs, the 
spectra have been shifted in flux in steps of half 
a decade. Different air shower components, in- 
cluding electrons, muons, hadrons, Cerenkov, 
and fluorescence photons, have been used to de- 
rive the spectra, see Table ||, and the results are 
compatible with the average all-particle spec- 
trum. However, for some experiments system- 
atic discrepancies to the all-particle spectrum 
can be recognized. 

As the final result of the present investigations 
the all-particle energy spectrum derived accord- 
ing to the poly-gonato model with a rigidity de- 
pendent cut-off and a common A7 is given in 
Table || for reference. The differences between 
the calculated spectrum and the average mea- 
sured flux 
errors. 



see Figure lib — are given as 
The measurements start at 0.5 PeV in 
the figure, consequently, below this energy no 
errors are quoted. The mean relative deviation 
of the calculated spectrum amounts only to re- 
markable 3%. 

6 Cosmic— ray mass composi- 
tion 

An often-used quantity to characterize the 
cosmic-ray mass composition above 1 PeV is the 
mean logarithmic mass, defined as 



{In A) =J2nlnAi 



(7) 



elements with the non-linear dependence of 7^ 
on Z are important to describe the measured 



T{ being the relative fraction of nuclei of mass A(. 
Some experiments derive this observable under 
certain assumptions from their measurements. 
In the superposition model of air showers, the 
shower development of heavy nuclei with mass A 
and energy Eq is described by the sum of A pro- 
ton showers of energy E = Eq/A. The shower 
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Figure 16: Average depth of the shower maximum X max versus primary energy as obtained by different 
experiments, for references see Table ||. Simulations taken from |112| , 113] using CORSIKA/QGSJET (solid 



lines) and MOCCA (dashed lines) are given for primary protons and iron nuclei. 



maximum t penetrates into the atmosphere as 
t oc InE, hence, most air shower observables at 
ground level scale proportional to In A Several 
experiments publish observables, from which the 
mean logarithmic mass can be derived. 

6.1 Average depth of shower maxi- 
mum 

A typical observable is the average depth of the 
shower maximum X max in the atmosphere. Fig- 
ure 16 summarizes recent measurements using 



Cerenkov and fluorescence light observations. In 
addition, results from Haverah Park are pre- 
sented. This experiment measures the electro- 
magnetic component only, no Cerenkov or fluo- 
rescence light. X max is derived from the struc- 
ture of the arrival time distribution and the lat- 
eral distribution of secondary particles. Charac- 
teristics of the experiments are compiled in the 
first part of Table || 

The results show systematic differences of « 
30 g/cm 2 at 1 PeV increasing to ~ 65 g/cm 2 
close to 10 PeV. The strongest increase of X max 
as function of energy is reported by the DICE 
experiment. A recent reanalysis indicates, that 
part of this behaviour could be caused by sys- 



tematic effects, including electronic saturation 
in the photomultiplier tubes, which had not 
been taken into account in the original analy- 
sis iTTc 



Only the Fly's Eye, HiRes, and DICE exper- 
iments observe an image of the air shower and 
derive the average depth of the shower maxi- 
mum directly. All other experiments use non- 
imaging techniques, i.e. the average depth of 
the shower maximum is derived from properties 
of secondary particle distributions at ground. 
Hence, the results obtained depend on the mod- 
els used to describe the air shower develop- 
ment. The programs applied by different groups 
are listed in Table || Many groups use COR- 
SIKA with the interaction models QGSJET or 
VENUS, but AIRES, KNP, and MOCCA are 
also applied. A private model including mini- 
jets is used for the Mt. Lian Wang analysis. Ap- 
plication of different codes can explain parts of 
the discrepancies in X max . For example, Dick- 
inson et al. found that using different interac- 
tion models in the MOCCA code, viz. the origi- 
nal and the SIBYLL model, the systematic error 
amounts to AX max m 10 g/cm 2 |]102[|. 



The measurements are compared to predic- 



17 



Table 5: Experiments and their corresponding air shower components analyzed to derive the mean logarithmic 
mass: e: electromagnetic, fj,: muonic, h: hadronic component, C: Cerenkov, and F: fluorescence light. The 
energy threshold for muons and hadrons is given. The last column lists the model used to interpret the data. 



Experiment 



H (E 



th) 



C F model to derive X„ 



BLANCA 
CACTI 
DICE (7| 
Fly's Eye @ 
Haverah Park 
HEGRA (Airobic) ]T5 
HiRes/MIA g0(| 
Mt. Lian Wang fTM]] 
SPASE/ VULCAN M 
Yakutsk fl03] 



x x CORSIKA QGSJET 

x x MOCCA-92 

x 800 MeV x CORSIKA VENUS 
x KNP 

x AIRES/SIBYLL 

x x CORSIKA QGSJET 

800 MeV x CORSIKA QGSJET 

x x MC with minijet model 

x x MOCCA SIBYLL 

x x QGSJET 



Experiment 



th i 



h (E 



th) 



air shower model 



CASA-MIA pi 
Chacaltaya |105| 1 
EAS-TOP/MACRO p06| 
EAS-TOP p7fl 
HEGRA (CRT) p8[ 
KASCADE (hadrons/muons) g09| 
KASCADE (electrons/muons) @] 
KASCADE (neural network) |l| 



x 800 MeV MOCCA -> QGSJET 

x 5 TeV CORSIKA QGSJET 

x 1.3 TeV CORSIKA QGSJET 

x 1 GeV CORSIKA QGSJET 

x track angle CORSIKA VENUS 

230 MeV 50 GeV CORSIKA QGSJET 

x 230 MeV CORSIKA QGSJET 

x 230/2400 MeV 100 GeV CORSIKA QGSJET 



tions of the shower maximum for protons and 
iron nuclei from simulations by Heck et al., using 
CORSIKA/QGSJET (ll2|, and by Pryke, using 
the MOCCA code with the internal hadronic in- 
teraction model |l|]. In the MOCCA model 
the energy is transported deeper into the atmo- 
sphere and the shower maxima are shifted by 
about 50 g/cm 2 to larger X max values as com- 
pared to CORSIKA/QGSJET. This results in 
a heavier composition. If the mean logarith- 
mic mass is derived from X max according to the 
MOCCA simulations, the data above 10 PeV 
indicate a pure Fe composition. As mentioned 
above CORSIKA/QGSJET seems to describe 
the development of air showers best. Hence, in 
the following the QGSJET interpretation is used 
to derive (lnj4). 

Despite of the differences mentioned, a gen- 
eral trend is visible in the data. All results start 
with an intermediate composition. However, the 



measured X max increase faster as function of en- 
ergy than the model predictions, i.e. the mea- 
sured values approach the proton distribution 
up to about 4 PeV. Above this energy, a change 
in the shape of the measured X max curves can 
be recognized in Figure ^ and the data start to 
approximate the calculations for iron nuclei with 
a nearest approach at about 30 PeV. Finally, at 
higher energies the measurements move again in 
direction of the proton predictions. 

Knowing the average depth of the shower 
maximum for protons X'm ax and iron nuclei 
X^ ax from simulations, the mean logarithmic 
mass can be derived in the superposition model 
of air showers from the measured X™^x S using 



(In A) 



vmeas vp 



X Fe 



■In A 



Ft 



(8) 



The corresponding (In A) values, obtained from 



the results shown in Figure 16, are plotted in 
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Figure 17: Mean logarithmic mass vs. primary energy, a) Results from the average depth of the shower 
maximum X max using CORSIKA/QGSJET simulations, b) Results from measurements of distributions for 
electrons, muo ns, a nd hadrons at ground level. For references see Fable g. Results from the balloon experi- 
ments JACEE |l!4| and RUNJOB are given as well. Predictions according to the poly-gonato model are 
represented by the solid line. Fhe dashed lines are obtained by introducing an ad-hoc component of hydrogen 
only, see text. 



Figure 17 a versus the primary energy. One ob- 
serves a large scattering of the values. In the 
region 1 to 10 PeV the maximum variations 
amount to A (In A) ~ 2. 



6.2 Particle distributions at ground 
level 



Results on (In A) from experiments measur- 
ing the electromagnetic, muonic, and hadronic 
components at ground level are shown in Fig- 
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ure 17 b. The experiments and the observables 
used as well as the air shower model applied 
are listed in the second part of Table pi The 
number of electrons and the number of muons 
are used by CASA-MIA, EAS-TOP/MACRO, 
and KASCADE (electrons/muons) . EAS-TOP 
uses the number of muons at a distance of 
200 m to the shower core in addition to the 
number of electrons. The hadronic component 
is analyzed by the Chacaltaya and KASCADE 
groups. HEGRA (CRT) measures the longitudi- 
nal shower development using the average pro- 
duction height of muons, i.e. the angle of inci- 
dence for individual muon tracks. Most experi- 
ments use parametric methods but also classifi- 
cations on an event-by-event-basis are applied, 
e.g. in the KASCADE neural network analysis. 

Several groups study the systematic effect 
caused by different interaction models. If pos- 
sible, the CORSIKA QGSJET interpretation is 
utilized for the compilation. CASA-MIA uses 
MOCCA/SIBYLL simulations but in their anal- 
ysis the results have been modified in order to 
match results from QGSJET p4j| . 

In the KASCADE neural network analysis 
[ jSl| l several combinations of observables and 
their influence on (In A) have been investigated. 
The observables include the number of electrons, 
the number of muons for energy thresholds of 
230 MeV and 2.4 GeV, the number of hadrons 
with Eh > 100 GeV and their energy sum as 
well the spatial structure of high-energy muons. 
In the figure two extreme cases are shown with 
a difference of about 0.7 to 1.5 in (In A). The 
lower values are obtained by analyzing the num- 
ber of electrons and the number of high-energy 
muons, while the high values result from investi- 
gations of the number of low-energy muons, the 
hadronic energy sum, and the spatial structure 
of the high-energy muons. 

In another investigation of several hadronic 
observables by the KASCADE group the sys- 
tematic differences amount to A In .A ~ 0.4 
[109], compatible with an estimate of Swordy 
et al. claiming an uncertainty of a {in A) ~ 0.23 
[111]. This, most likely, represents the system- 



In both graphs of Figure 17 the mean loga- 
rithmic mass calculated according to the poly- 
gonato model with rigidity dependent cut-off is 
given by the solid line. Since more and more el- 
ements reach their cut-off energy Ez, the mean 
logarithmic mass increases with rising energy 
and would finally reach pure uranium. 



As can be seen in Figure 14 and has been dis- 
cussed in section 



above 100 PeV the sum 
spectrum for 1 < Z < 92 does not sufficiently 
describe the all-particle spectrum. For this rea- 
son, a new component is introduced ad hoc in 
order to fill the difference between the sum spec- 
trum (solid graph) and the all-particle spectrum 



represented by the dashed line in Figure 14. Sev- 



eral theories predict an extragalactic component 
consisting of hydrogen and helium nuclei only, 
see e.g. [10|. In fact, the AKENO experiment 
[ ]94| found that the mass composition becomes 
heavier with energy up to values of (In A) ~ 3 
at iV e ~ 10 7 ' 5 . For larger showers the composi- 
tion gets lighter again, and in analyses with two 
or three mass groups the proton component re- 
covers. The conversion of N e to primary energy 
had not been done, but cautious estimations re- 
late lgiV e = 7.5 to about 300 PeV. 

Anticipating a pure proton composition yields 
a lower limit for the (In A) values independent 
of any model assumptions. Hence, only protons 
are assumed for the ad-hoc component. The 
mean logarithmic mass obtained is presented in 



Figure 17 as dashed lines. 



atic uncertainties introduced by inconsistencies 
in the hadronic interaction models and M.C. 
programs used to interpret the data. 



The results from experiments measuring elec- 
trons, muons and hadrons at ground level are in 
reasonable agreement with direct measurements 
in the region of overlap, namely below 1 PeV. Up 
to about 3 PeV they show a small trend of an 
increasing (In A). For higher energies the mea- 
sured (In A) values indicate a more pronounced 
increase. Around 30 PeV values of (In A) « 3.3, 
corresponding to silicon, are obtained. The pre- 
diction of the poly-gonato model yields the same 
shape of (In A) as function of energy. On the 
other hand, systematic differences between the 
results derived from the X max observations and 
the model predictions are apparent. 
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6.3 X max versus ground level observ- 
ables 

In order to study the difference in shape and 
absolute value of (In A) in more detail, the 
weighted average of (In A) is calculated for both, 
ground observables and X max measurements, 
taking into account the given errors for the 
individual measurements. For the calculation 
of (In A) from X max simulations with COR- 
SIKA/QGSJET fTT^] and the MOCCA code 
[|113| are used. Both results are shown in Fig- 
ure [0]a together with the corresponding aver- 
age values for ground observables in Figure |l£|b. 
The contour lines represent the r.m.s. values of 
the data. 

The parameters of the poly-gonato model 
have been obtained by fits to the all-particle 
energy spectrum. The mean logarithmic mass 
is an independent observable, which allows an 
additional check of the model. The values pre- 
dicted by the model are compared to the data 



in Figure 15. To estimate the systematic uncer- 
tainties introduced by the extrapolation of 7^ 
for heavy nuclei, the mean logarithmic mass is 
shown for the non-linear (C > 1) and the linear 
(C = 1) extrapolation by the solid and dash dot- 
ted lines, respectively. In addition, (In A) values 
including an ad-hoc component of protons only 
are represented for both cases by the dashed and 
dotted lines, respectively. The differences be- 
tween the two extrapolations for 7^ are relative 
small, e.g. the maximum (In .A) values for the 
two cases including the ad-hoc component dif- 
fer only by A (In A) « 0.2. 

The variations of (In A) derived from ground 
observables are in the same order as the uncer- 
tainties due to the observables used. The in- 
vestigations by the K AS CADE group — as dis- 
cussed above — suggest a systematic error of 
about 0.5 to 0.7 for the determination of (In A) 
compatible with the average variation of ~ 0.5 
in Figure [l8|b. The results of the poly-gonato 
model are in reasonable agreement with the 
measurements, absolute values and the shape 
are almost the same. Even the maximum of 
(In A) around 70 PeV is reflected by the mea- 
surements. Below 0.5 PeV only a few indirect 
measurements contribute, see Figure 17, and 
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Fi gure 18: Average mean logarithmic mass (lnA) 
vs. primary energy. Predictions of the poly-gonato 
model are given without and with the ad-hoc com- 
ponent as well as for the non-linear (C > 1) and 
linear (C = 1) extrapolation of 72. They are com- 
pared to averaged values for experiments: a) ob- 
serving X max using two different interaction models 
for the interpretation, b) measuring electrons, muons 
and hadrons at ground level. The contour lines rep- 
resent the r.m.s. values of the data. 



poly-gonato model are not to be taken seriously. 
Above 20 PeV the model yields a heavier com- 
position as compared to the mean value of the 
measurements. It can not be excluded, that in 
some analysis procedures the mass values had 
been restricted to a range limited by the ex- 
treme values pure protons and pure iron nuclei. 
Therefore, the numbers at high energies could 
be biased towards a lighter composition. 

The mean logarithmic mass derived from 
shower maximum observations shows for both 
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simulation codes below the knee energy a 
slightly falling (In A) . Beyond this energy (In A) 
increases, but less pronounced as compared to 
the ground observations. The (In A) values 
reach their maximum at about 50 PeV. This 
energy is in agreement with the phenomeno- 
logical model, where the maximum is caused 
by the introduction of the ad-hoc component. 
As can be inferred from Figure 18 a, the re- 
sults of the X max measurements are not com- 
patible with the poly-gonato model. Especially 
below the icnee they do not match the results 
of direct observations, nor their extrapolations 
to higher energies. The decreasing (In A) val- 
ues below the knee are not seen by direct mea- 
surements. This suggests a systematic difference 
between ground observables and X max measure- 
ments which could have its origin in the air 
shower simulations or data analysis procedures 
applied. The effect is independent of the sim- 
ulation program used. Both, CORSIKA and 
MOCCA, yield different absolute values, but the 
same inclined "S" shape is obtained. 

6.4 Relation for cut— off energy 

The mean logarithmic mass calculated accord- 
ing to the poly-gonato model with mass depen- 
dent cut-off shows a similar shape but yields 
larger values as compared to a rigidity depen- 
dent cut-off. The results for a common 7 C using 
the non-linear extrapolation for jz are plotted 
in Figure ^ versus the particle energy and com- 
pared to the average results of the experiments. 
While the (In A) values for the rigidity depen- 
dent cut-off are very close to the average mea- 



sured results, see Figure |18|b, the values for a 
mass dependent cut-off above 1 PeV are about 
0.3 larger than the average data. Thus, a rigidity 
dependent cut-off is favoured against the mass 
dependent approach. The hypothesis of a con- 
stant cut-off energy yields an almost constant 
mean logarithmic mass, as shown in Figure 
and is not compatible with the measurements. 

Combining both, the arguments discussed 
above in connection with the shape of the all- 
particle energy spectrum and the findings con- 
cerning the mean logarithmic mass, the poly- 
gonato model with rigidity dependent cut-off, 
the non-linear extrapolation of jz, an d a com- 
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Figure 19: Average mean logarithmic mass (\nA) 
vs. primary energy. Predictions of the poly-gonato 
model with mass dependent cut-off (solid line) and 
with constant cut-off (dashed line) are compared 
to averaged experimental values. The experimental 
data in a) and b) are the same as in Figure 18. 



mon A7 is the favoured solution. The (ln^4) 
values versus energy for this model are summa- 
rized in Table |9] for reference. 

7 Compatibility with results 
from indirect measurements 

The parameters of the poly-gonato model have 
been determined by extrapolating the directly 
measured energy spectra to high energies and 
a fit to the all-particle spectrum of indirect 
mesurements in section ||. As a first indepen- 
dent test, the mean logarithmic mass calculated 
with the poly-gonato model has been compared 
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to experimental results in section ||. In the fol- 
lowing, additional checks of consistency are pre- 
sented. Two results from air shower measure- 
ments are compared to the predictions of the 
model. Namely the spectra for mass groups ob- 
tained by the K AS CADE experiment |7^] and 
Xmax distributions published by the Fly's Eye 



group [115 1 to test the predictions in the energy 
regions of the knee and the cut-off for the ultra- 
heavy elements, respectively. 



7.1 Energy spectra for groups of ele- 
ments 

Recently, individual spectra for groups of el- 
ements in the energy range from 10 6 GeV 
to 10 8 GeV have been reconstructed from air 
shower measurements by the K AS CADE group 
Using shower size spectra for the elec- 
tromagnetic and muonic component for three 
different zenith angle bins the primary energy 
spectra for elemental groups have been obtained 
by an unfolding procedure. Four groups, rep- 
resented by hydrogen, helium, CNO, and iron 
have been chosen. These groups are identified 
with the following groups in the poly-gonato 
model: Z = 1, 2 < Z < 5, 6 < Z < 24, and 
25 < Z < 92. Renormalizing the KASCADE 
energy scale (according to Table || by -7%) the 
spectra shown in Figure |20] are obtained. 



> 
a> 

CD ? 
,-10 3 

'w 



"^B.g.Pap 



1 

2-5 
6-24 
25-92 



w 10 

o 
LU 

o 

T3 
X 
ZJ 




« 



O p 

- □ He 

z A CNO 

: Fe 

• all 

II I L 



10 



10 10 
Energy E [GeV] 



Figure 20: Energy spectra for four elemental groups 
and the all-particle spectrum as measured by the 
KASCADE experiment [[79) (symbols) compared to 
the spectra of the poly-gonato model (lines). 



The measured all-particle spectrum agrees 
well with the calculated flux also shown in the 
figure. The shapes of the experimental spec- 
tra for protons and the helium group are com- 
patible with the spectra calculated according to 
the poly-gonato model. It is remarkable that 
the cut-off behaviour of both components agrees 
well with the calculated spectra, thus indepen- 
dently confirming the cut-off parameters ob- 
tained in section |5[ The measured spectra for 
the CNO group and the iron group also agree 
with the calculations, with somewhat larger ex- 
perimental error bars at high energies. In the 
energy region around 10 6 GeV the KASCADE 
fluxes for protons and helium are slightly larger 
than in the poly-gonato model and the values 
for the iron group are somewhat below the cal- 
culated flux. On the whole, in an overall view 
the spectra of the poly-gonato model agree well 
with the unfolding analysis. It may be pointed 
out that both analyses are completely indepen- 
dent and they yield very similar results. 



7.2 



X max — distributions 



The energy region between 10 8 GeV and 
10 9 GeV where the heaviest elements reach 
their cut-off energies and the ad-hoc compo- 
nent starts to dominate the all-particle spec- 
trum is an interesting range for further indepen- 
dent tests of the predictions of the model. 

Assuming protons only for the ad-hoc compo- 
nent, relative abundances for elemental groups 
are obtained as compiled in Table || for energies 
E = 3 • 10 8 GeV, 5 • 10 8 GeV, and 10 9 GeV. 
In the energy range mentioned the fraction of 
heavy and ultra-heavy elements is expected to 
decrease from 62% to 35%. On the other hand 
the fraction of the ad-hoc component steadily 
increases, accounting for the rise of the light el- 
ements (protons). 

Adopting the mass composition listed in Ta- 
ble ^, simulations have peen performed in or- 
der to study the expected distributions of the 
average depth of the shower maximum X max . 
The air shower simulation program CORSIKA 
(6.0120) with the interaction model QGSJET 
has been used to calculate X max for primary 
protons and iron induced showers for the ener- 
gies listed in the table. For each energy and 
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Table 6: Relative abundances for groups of nuclei 
with charge number Z for different primary energies 
Eq according to the poly-gonato model. The ad-hoc 
component is assumed to be protons only. 



0.01 







Energy Eo/GeV 




Z 


3 • 10 s 5 • 10 8 10 9 


light 


1-12 


38% 52% 65% 


heavy 


13-29 


24% 12% 8% 


ultra heavy 


30-92 


38% 36% 27% 



species 150 showers have been generated. The 
CORSIKA code allows primary nuclei with mass 
numbers A < 56. The X max distributions for 
ultra-heavy nuclei have been estimated using 
the assumption, X max oc In A, based on the sim- 
ulated values for protons and iron nuclei. 

In section |0] it has been discussed, that there 
are discrepancies in the mean logarithmic mass 
calculated from electron, muon, and hadron dis- 
tributions on one hand and X max values on the 
other hand. Therefore, it can not be expected, 
that the CORSIKA predictions match the mea- 
sured values. The distances between the lines 
in Figure ^ suggest, that the differences be- 
tween the average X max values for protons and 
iron nuclei are almost independent from the in- 
teraction model. If a constant offset is added 
to the simulated values, only the absolute val- 
ues of the average depths of the shower maxi- 
mum are shifted but the fluctuations stay the 
same. Hence, to reconcile the simulations with 
the measured data, an offset AX max = 30 g/cm 2 
has been added to the simulated values. The re- 
sulting distributions are compared in Figure ^ 
to measurements of the Fly's Eye group [115]. 

As can be inferred from the figure, with this 
offset the simulated and measured distributions 
are compatible with each other. For the low- 
est energy bin the assumption of two groups 
only for elements up to iron is only a crude ap- 
proximation of the real mass composition and 
the simulation of intermediate elements like the 
CNO group would improve the agreement be- 
tween the simulated and measured distributions. 
The structures on the right hand side of the pro- 
ton distributions are due to fluctuations in the 
simulated X max values. 

The measured showers seem to penetrate 
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Figure 21: Distributions for the average depth of the 
shower maximum X max for different energy ranges: 3 
to 5 -10 s GeV (top), 5-10 8 GeVto 10 9 GeV (middle), 
and above 10 9 GeV (bottom) Th e measured values 
of the Fly's Eye experiment [115| are compared to 
simulated values for protons, iron and ultra-heavy 
nuclei, see text. 



deeper into the atmosphere than predicted by 
the CORSIKA calculations but the width of the 
measured distributions, i.e. the X max fluctua- 
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tions are well compatible with the mass compo- 
sition predicted by the poly-gonato model. 



8 Conclusion 

An all-particle primary energy spectrum has 
been combined from many air shower experi- 
ments by properly adjusting the individual en- 
ergy scales to meet the flux values of direct mea- 
surements. The resulting power law spectrum 
exhibits two changes in the spectral slope at 
4.5 PeV and about 0.4 EeV. 

A phenomenological model, named poly- 
gonato model, is adopted in which the energy 
spectra of individual nuclei have a power law be- 
haviour with a cut-off at a specific energy Ez- 
The icnee is explained as the subsequent cut-offs 
of the individual elements of the galactic com- 
ponent, starting with protons. The second knee 
seems to indicate the end of the stable elements 
of the galactic component. 

It turns out, that the contribution of ultra- 
heavy elements with Z > 28 is not negligible in 
the 100 PeV regime. Their flux values have been 
calculated, taking into consideration the modu- 
lation in the heliosphere and power law spectra, 
for which the spectral indices have been esti- 
mated in two independent ways. 

Several assumptions about the cut-off be- 
haviour for the individual elements have been 
checked against the experimental data. It turns 
out that the hypothesis the icnee energy scales 
with the charge of the individual nuclei, Ez = 
Z ■ E p , describes the data best. 

Comparing several methods of extracting a 
mean logarithmic mass from the data indicates 
a systematic error of A (In A) ss 0.5 for the de- 
termination of (In A) . Taking this error into ac- 
count, the mass composition calculated with the 
poly-gonato model is in agreement with results 
from air shower experiments measuring the elec- 
tromagnetic, muonic and hadronic components 
at ground level. But the mass composition dis- 
agrees with results from experiments measuring 
the average depth of the shower maximum with 
Cerenkov and fluorescence detectors. 
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Table 7: Absolute flux [(m 2 sr s TeV) 1 ] at Eq = 1 TeV/nucleus and spectral index 72 of cosmic-ray 
elements. 
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7.79 


10" 


-7 


2.46 


73 4 


Ta 


2.16 


10" 


-7 


2.16 


12 3 


Mg 


8.01 


10- 


-3 


2.64 


43 4 


Tc 


6.98 


10" 


-8 


2.46 


74 4 


W 


4.16 


10- 


-7 


2.15 


13 3 


Al 


1.15 


10" 


-3 


2.66 


444 


Ru 


3.01 


10- 


-7 


2.45 


75 4 


Re 


3.35 


10- 


-7 


2.13 


143 


Si 


7.96 


10" 


-3 


2.75 


45 4 


Rh 


3.77 


10- 


-7 


2.44 


76 4 


Os 


6.42 


10- 


-7 


2.12 


15 3 


P 


2.70 


10" 


-4 


2.69 


46 4 


Pd 


5.10 


10- 


-7 


2.43 


77 4 


Ir 


6.63 


10- 


-7 


2.11 


16 3 


S 


2.29 


10" 


-3 


2.55 


474 


Ag 


4.54 


10- 


-7 


2.42 


78 4 


Pt 


1.03 


10- 


-6 


2.10 


17 3 


CI 


2.94 


10" 


-4 


2.68 


48 4 


Cd 


6.30 


10- 


-7 


2.41 


79 4 


Au 


7.70 


10- 


-7 


2.09 


18 3 


Ar 


8.36 


10" 


-4 


2.64 


49 4 


In 


1.61 


10" 


-7 


2.40 


80 4 


Hg 


7.43 


10" 


-7 


2.08 


19 3 


K 


5.36 


10- 


-4 


2.65 


50 4 


Sn 


7.15 


10" 


-7 


2.39 


81 4 


Ti 


4.28 


10- 


-7 


2.06 


20 3 


Ca 


1.47 


10" 


-3 


2.70 


51 4 


Sb 


2.03 


10" 


-7 


2.38 


82 4 


Pb 


8.06 


10" 


-7 


2.05 


21 3 


Sc 


3.04 


10- 


-4 


2.64 


52 4 


Tc 


9.10 


10" 


-7 


2.37 


83 4 


Bi 


3.25 


10- 


-7 


2.04 


22 3 


Ti 


1.14 


10" 


-3 


2.61 


53 4 


I 


1.34 


10- 


-7 


2.37 


84 4 


Po 


3.99 


10" 


-7 


2.03 


23 3 


V 


6.31 


10- 


-4 


2.63 


54 4 


Xe 


5.74 


10- 


-7 


2.36 


85 4 


At 


4.08 


10- 


-8 


2.02 


24 3 


Cr 


1.36 


10" 


-3 


2.67 


55 4 


Cs 


2.79 


10- 


-7 


2.35 


86 4 


Rn 


1.74 


10- 


-7 


2.00 


25 3 


Mn 


1.35 


10" 


-3 


2.46 


56 4 


Ba 


1.23 


10- 


-6 


2.34 


87 4 


Fr 


1.78 


10- 


-8 


1.99 


26 2 


Fe 


2.04 


10" 


-2 


2.59 


57 4 


La 


1.23 


10- 


-7 


2.33 


88 4 


Ra 


7.54 


10" 


-8 


1.98 


27 3 


Co 


7.51 


10- 


-5 


2.72 


58 4 


Ce 


5.10 


10- 


-7 


2.32 


89 4 


Ac 


1.97 


10- 


-8 


1.97 


28 3 


Ni 


9.96 


10" 


-4 


2.51 


59 4 


Pr 


9.52 


10" 


-8 


2.31 


90 4 


Th 


8.87 


10- 


-8 


1.96 


29 4 


Cu 


2.18 


10" 


-5 


2.57 


60 4 


Nd 


4.05 


10- 


-7 


2.30 


91 4 


Pa 


1.71 


10- 


-8 


1.94 


30 4 


Zn 


1.66 


10" 


-5 


2.56 


61 4 


Pm 


8.30 


10- 


-8 


2.29 


92 4 


U 


3.54 


10" 


-7 


1.93 


31 4 


Ga 


2.75 


10- 


-6 


2.55 


62 4 


Sm 


3.68 


10" 


-7 


2.28 















2 This work, see Figures [2] to and text. 
3 From Wiebel-Sooth et al. fl6|| . 

4 This work, extrapolation for ultra-heavy elements, see text. 
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Table 8: All-particle cosmic-ray energy spectrum. Normalized average flux from air shower experiments $ m 
and flux <!> c according to the poly-gonato model. The errors are r.m.s. values of the data for $ m and the 
difference between $ m and $ c for $ c . {lgE [GeV] and $(£ ) ■ Efr 5 [GeV 1 ' 5 /™ 2 sr s].) 









$ c (£o) 


P 2.5 
• Eq 






• E^ 


<S>c(E ) 




4.0 






4974 




6.0 


2417± 


207 


2395± 


22 


4.1 






4796 




6.1 


2270± 


212 


2292± 


21 


4.2 






4624 




6.2 


2164± 


177 


2182± 


19 


4.3 






4459 




6.3 


2051± 


191 


2065± 


15 


4.4 






4300 




6.4 


1951± 


203 


1939± 


12 


4.5 






4146 




6.5 


1843± 


230 


1802± 


40 


4.6 






3999 




6.6 


1680± 


236 


1658± 


22 


4.7 


3911± 


194 


3857± 


54 


6.7 


1519± 


245 


1509± 


10 


4.8 


3841± 


201 


3720± 


121 


6.8 


1368± 


250 


1362± 


7 


4.9 


3711± 


197 


3589± 


122 


6.9 


1229± 


213 


1220± 


9 


5.0 


3579± 


194 


3463± 


116 


7.0 


1076± 


204 


1090± 


14 


5.1 


3427± 


192 


3341± 


86 


7.1 


lOllzb 


229 


974± 


38 


5.2 


3345± 


193 


3224± 


121 


7.2 


914± 


262 


871± 


43 


5.3 


3218± 


191 


3111± 


108 


7.3 


815± 


333 


782± 


33 


5.4 


3086± 


201 


3002± 


85 


7.4 


666± 


187 


702± 


36 


5.5 


2917± 


251 


2896± 


21 


7.5 


586± 


115 


629± 


43 


5.6 


2772± 


252 


2794± 


22 


7.6 


514± 


97 


559± 


46 


5.7 


2704± 


237 


2693± 


11 


7.7 


461± 


100 


493± 


33 


5.8 


2645± 


235 


2594± 


50 


7.8 


426± 


95 


429± 


4 


5.9 


2526± 


217 


2495± 


31 


7.9 


353± 


75 


369± 


16 












8.0 


323± 


54 


314± 


9 



Table 9: Mean logarithmic mass according to the poly-gonato model (lg£b [GeV]). 



Ig^o 


(In A) 


lg£o 


(In A) 


lg£o 


(In A) 


lg£o 


(In A) 


4.0 


1.66 


4.1 


1.68 


4.2 


1.69 


4.3 


1.70 


4.4 


1.72 


4.5 


1.73 


4.6 


1.74 


4.7 


1.75 


4.8 


1.77 


4.9 


1.78 


5.0 


1.79 


5.1 


1.81 


5.2 


1.82 


5.3 


1.83 


5.4 


1.85 


5.5 


1.86 


5.6 


1.88 


5.7 


1.90 


5.8 


1.92 


5.9 


1.94 


6.0 


1.97 


6.1 


2.00 


6.2 


2.04 


6.3 


2.09 


6.4 


2.16 


6.5 


2.24 


6.6 


2.35 


6.7 


2.47 


6.8 


2.61 


6.9 


2.77 


7.0 


2.93 


7.1 


3.09 


7.2 


3.24 


7.3 


3.38 


7.4 


3.51 


7.5 


3.62 


7.6 


3.72 5 -3.73 6 


7.7 


3.80-3.83 


7.8 


3.87-3.92 


7.9 


3.92-4.02 


8.0 


3.93-4.12 


8.1 


3.88-4.24 


8.2 


3.73-4.36 


8.3 


3.50-4.50 



5 1 < Z < 92 + ad-hoc component, protons only. 
6 1 < Z < 92 only. 
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